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Abstract 

A low order diagrammatic study of the dimension dependent Su-Schrieffer- 
Heeger model Hamiltonian in the weak electron-phonon coupling regime is 
presented. Exact computation of both the charge carrier effective mass and 
the electron spectral function shows that electrons are good quasiparticles in 
the antiadiabatic limit but new features emerge in the adiabatic and inter- 
mediate regime, where the phonons and the electrons compare on the energy 
scale. Here we find: i) a sizeable mass enhancement over the bare band value, 
ii) the appearance of many transition peaks in the band bottom spectral func- 
tion together with a growing loss of spectral weight at larger e-ph couplings. 
The onset of a polaronic state is favoured in two dimensions. 
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l.INTRODUCTION 



A sizeable electron-phonon interaction can induce a local deformation in the lattice 
structure accompanied by the formation of a quasiparticle with multiphononic character, 
the polaron |IJ]. A wide literature has been produced on the subject during the last decades 
and particular emphasis has been recently laid on the polaronic properties of high 
T c superconductors ||. As the spatial extension of the lattice deformation can vary, the 
concepts of large and small polaron have been introduced: the transition between a large 
and a small polaron state is driven by the strength of the electron-phonon coupling [pH|-rr3|] 
and monitored through the behavior of ground state properties such as the polaron energy 
band and the effective mass |T^|T6|. This transition [ T7j , leading to a self trapped state |18 



at strong couplings, is accompanied by a sizeable enhancement of the effective mass whose 
value may change considerably according to the degree of adiabaticity and the peculiarities 
of the lattice structure |T9|j20[| . Theoretical investigations usually start from the Holstein 
molecular crystal model jH]] which assumes a momentum independent coupling of electrons 
to dispersive optical phonons. The coupling to acoustical phonons, although possible in 
principle ||20|| , would lead to huge mass renormalizations ||22|| . In fact, also the Holstein 
optical polaron masses are very heavy (at least larger than 10 3 times the bare band mass) 
in the self-trapped state but masses of order 10 times the bare band mass are possible in 
the presence of high energy phonon spectra p3fl . As a remarkable feature of the Holstein 
model the polaron mass turns out to be essentially dimension independent for any value of 



the adiabaticity parameter [^4 



In some systems however [pq -|27| the e-ph interaction modifies the electron hopping 
matrix elements thus leading to a momentum dependent coupling function. In this case an 
appropriate theoretical framework is offered by the Su-Schrieffer-Heeger(SSH) model Hamil- 
tonian p8| originally proposed to explain the conducting properties of quasi one dimensional 
polymers as polyacetylene ||29|| . In these systems the CH monomers form chains of alternat- 
ing double and single bonds leading to two regions with different structural patterns having 
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the same energy. This twofold degenerate ground state sustains a nonlinear localized exci- 
tation, a domain wall separating the two regions. Being the wall thickness much larger than 



the lattice constant in the weak coupling regime a continuum version |3(| of the SSH model 
has been formulated thus leading to an analytical description of the solitonic excitation in 
the dimerized system The continuum model admits also polaron-like solutions p^-|3~l| 
and a periodic array of solitons solution, the soliton lattice |]35| , |36| . 

In general, the SSH Hamiltonian provides an alternative (to the Holstein model) tool to 
analyse the physics of polaron formation as tuned by the strength of the e-ph coupling. In 
two dimensions, the SSH model has also been studied as a particular case of the 2D half 



filled Hubbard model (with zero on site repulsion) [37| soon after the discovery of high T c 



superconductivity. In particular, a recent investigation |38| has revealed the complexity of 
the square lattice structure described by a SSH type Hamiltonian pointing out that the 
opening of the gap at the Fermi level (due to the Peierls distortion) involves many lattice 
modes with wave numbers parallel to the nesting vector. The static and dynamical polaronic 
properties have been also analysed in the same model and estimates of the effective mass 
have suggested that 2D polarons are heavier than ID polarons. 

This paper deals with the two dimensional electron-lattice system treating the SSH tight 
binding Hamiltonian by a weak coupling perturbative method. This approach, although not 
adequate to capture the full multiphononic nature of the polaronic quasiparticle, still can 
provide useful informations |40| regarding the onset of polaron formation in some portions of 
parameter space. Here we look first at the mass renormalization in one and two dimensions 
exploring a wide range of values for the adiabatic parameter and, successively, we compute 
the electronic spectral function to detect whether and to which extent bare electrons behave 
as good quasiparticles. The Section 2 outlines the SSH model and contains the results of 
this study while some conclusions are drawn in Section 3. 
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2. MODEL AND RESULTS 



In real space the SSH Hamiltonian reads 



H= E ^r + s(/ r f /r + s + fLfr) +E(|^+E f K " «r+s) 2 ) 

«/r,r+s = -^[^ + «K - M r+S )] (0.1) 

where the double summation over r and s runs over first neighbors lattice sites. J is 
the nearest neighbors hopping integral and isotropic conditions are assumed in the square 
lattice, a is the electron-phonon coupling, u r is the dimerization coordinate which specifies 
the displacement of the r— lattice site from the equilibrium position, p r is the momentum 
operator conjugate to u r , M is the ion (ionic group) mass, K is the effective spring constant, 
fl and f r create and destroy electrons on the r— site. Let's expand the lattice displacement 
and its conjugate momentum in terms of the phonon creation and annihilation operators 6 q 
and 6 q and Fourier transform the electron operators 

u r = E / 1 0-q + K) exp(iq • r) 



Pr = i E y "^r( & -q ~ 6 i) ex P(^ ' r ) 

fr= ^E eX P(^ k - r )/k (0.2) 



N k 

in order to obtain the SSH Hamiltonian in momentum space: 
H = if o + H int 

H0= E^/k + E^q fo q 
k q 

H int = E»( k + <l» k )( 6t -q + 6q)/I+q/k 
k,q 

£ k = — Jcos(k • a) 

2 if . 2 /q • a\ 
uoi = 4 — sm 2 - — 
q M \ 2 J 

#(k + q,k) = , (sin((k + q) -a) -sink-a) (0.3) 

y/2MNuq 



where N is the total number of lattice sites, a = |a| is the lattice constant and the 
2D reduced Brillouin zone is spanned by the vectors hi = (2n/a,0) and b 2 = (0,27r/a). 
Being q = (q x ,q y ), the numerical integration constraint is defined by q x G [0, 27r/a], q y G 
[0,2ir/a — q x ]. In ID, the phonon dispersion relation is defined in the range q G [0, 7r] and, 
in the reduced Brillouin zone, the spectrum displays both an acoustic and an optical branch 
f4lfl . The model contains three free parameters: the hopping integral J, the zone boundary 



frequency u n = 2yK/M which coincides with the zone center optical frequency in the 
reduced zone scheme, the coupling constant a 2 /4K. 

The full electron propagator in the Matsubara Green's functions formalism is defined as: 

G(k,r) = - £(-l) n / dr 1 ...drJT r Ur)H int {T l ) ■ -^(r n )^(0)) n (0.4) 
n=o Jo X /0 

where (3 is the inverse temperature, T T is the time ordering operator, < ... >o indicates 

that thermodynamic averages are taken with respect to the unperturbed Hamiltonian and 

only different connected diagrams contribute to any order n. I have computed exactly 

the self-energy terms due to one phonon (n = 2 in eq.(4)) and two phonons (n = 4 in 

eq.(4)) scattering processes which determine the renormalized electron mass m e ff through 

the relations: 

m ef f __ 1 - d-ReE k (e)/de| k= o ; t= -j 
m 1 + 9i?eS k (e)/9£: k | k= o; £ =-j 

where, J ReS k (e) = ReT^\e) + ReT^ a \e) + ReT^ b \e) + ii!eE k 2c) (e) is the frequency 
dependent real part of the retarded self-energy. There are three contributions due to different 
connected two-phonons diagrams p2 |. Their effect is however confined to the intermediate 



regime in which u n is comparable to the electronic energy J. We set J = O.leV with the 
caveat that electron-electron correlations (weak in conducting polymers with wide 7r-electron 
bands and not taken into account by the SSH model) may become relevant in narrow band 
systems. This value is lower than those usually taken for the SSH adiabatic model but 
it allows us to discuss also a broad range of (anti) adiabatic parameters with reasonable 
choices of phonon energies. In the intermediate regime, the two phonons diagrams enhance 



the effective mass by ~ 15% with respect to the one phonon result. Instead, in the fully 
adiabatic and antiadiabatic regimes the two phonons contributions (evaluated at the band 
bottom) are negligible. Hereafter, the displayed results depend on the very one phonon self- 
energy E k ^(ie m ) term (e m = (2m + l)ir/ j3 with m integer number) whose finite temperatures 
analytic expression is given by: 



4%™)= -£s 2 (k,k-q) 



n B {.^a) + n F (-£ k _ q ) n B (^ q ) + n F (£ k _ q ) 



(0.6) 



Itm — ^k-q — ~ ^k-q + ^q 

riB and rip are the Bose and Fermi occupation factors respectively. 

Figure 1 shows, both in one and two dimensions, a sizeable mass enhancement in the 
intermediate regime with a pronounced spike at cu n ~ y/2J '. In 2D the effective mass is larger 
than in ID. The onset of a mass renormalization starting at cu^ ~ J/2 and, more evidently, 
at u n ~ J (together with the increased relevance of multiphonons contributions) signals that 
polaron formation is expected in this regime while no mass enhancement is obtained in the 
adiabatic and antiadiabatic limits. The same trend is observed both in ID and 2D. Let's 
analyse in detail the origin of the divergent-like mass behavior. The main contribution to 
both self-energy partial derivatives in eq.(5) comes from the lattice mode vectors connecting 
two electronic states such that J + e q — cj q ~ hence, from q- vectors satisfying the relation 

JV2 



\Jl — cos(q • a) 



Although for any value J < lo % a set of singular q- vectors does exist, their divergent contri- 
butions to the numerator and denominator in eq. (5) generally cancel out and no substantial 
effect is seen on the effective mass. Only in the case u n ~ \[2 J something special occurs 
due to scattering by phonons at the points such that \q x + q y \ = tt/2. Infact (9i?eSk(e) / de^ 
contains as modulation factor a cos(q • a) term which is peculiar of the band structure and 
responsible for the van Hove singularities in the density of states. The vanishing of this term 
at the same wave vectors which allow for energy conservation finally results in the abrupt 
increase of the effective mass observed in Fig.l. In the square lattice the divergent-like 
behavior is more evident since many points (only one point in ID) fulfill the simultaneous 
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occurence of the two singular effects. It should be emphasized that this phenomenon, rather 
than being a general feature of two dimensional systems, is related to the peculiarities of 
the electron band in the square lattice. 

Let's look now at the spectral function defined by A(k, e) = — 2/mG re j(k, e) to get more 
insight K3J into the suggestions proposed by the effective mass behavior. In terms of the 



retarded self-energy, obtained by eq.(6) through analytic continuation e m — > e + i8, A(k, e) 
reads: 

A(k, e) = 2n5(e - e k - ifeE k (e)) + R2/mS k (e) _ (Q ?) 

(e - e k - ReEJe)) + ImEJe) 



The first addendum in eq. (7) contributes when JmSk(e) = 0. The band bottom (|k| = 0) 
spectral function has been computed in a number of representative cases and numerical 
convergence has been achieved by summing the self-energy term (eq.(6)) over 6000 q— points 
in the ID Brillouin zone and 90000 q— points in the reduced 2D Brillouin zone. The 
complexity of the numerical work is mainly related to the search of the zeros of the 5- 
function argument in eq.(7) and in the im£fc=o(e). In ID, we have used the representation 

5(q - q°) 



sum 



\df/dq\ q=q o 



and obtained convergence by summing over 10 points of the e— axis. In 2D, the sum over 
q vectors can be handled as follows: 



dq y 



oo 

'""V 



N a N xa 27T 

with V being the cell volume and, at a fixed q x , the 5-function transforms as 

XUf \i S(q y -q°) 
S[f(q x ,q y )] = rfTTTTT , 

\ d j/ d( lyUy=q y 

At any e (we take 80000 points) the program searches the (q x , q y ) points which allow for 
energy conservation and the JmSfc =0 (e) is normalized over the total number of these pairs. 
The sum rule 

de 



/oo w f 
-A(\k\ = 0, 6) = 1 
-oo Z7T 



has to be numerically fulfilled in principle. Here however we are approximating the total self- 
energy by the one-phonon term (eq.(6)) which linearly depends on the free parameter a 2 /4K. 
As the e-ph coupling grows multiphonons terms become more relevant and our approximation 
becomes less accurate. Accordingly, deviations from the sum rule are expected as a measure 
of the loss of spectral weight associated with higher order self-energy effects. In this regard, 
the sum rule numerical analysis permits to define the range of a 2 /4K values within which 
the one phonon approximation is reliable. For any choice of input parameters we are able 
to estimate the intrinsic error of our physical model. 

Let's start the discussion looking at the spectral function in the adiabatic regime (Figures 
2). In one dimension and at very weak coupling (Fig. 2 (a)) the sum rule is satisfied and the 
spectral weight spreads mainly in a few peaks around the highest one located at the energy 
e = — lOOmeV with a significant tail up to energy levels of order e ~ — 80meV. At larger 
couplings (Fig.2(c)), where the dimensionless effective coupling is a 2 / (AKJ) = 0.1, the tail is 
appreciable up to e ~ — 50meV and the loss of spectral weight is ~ 20%, being /f^ deA(k = 
0, e)/27r = 0.8. A well defined transition is however still present at e = — 108.2meV. In 
two dimensions, a few well resolved transitions around e = — lOOmeV show up at very weak 
couplings (Fig. 2(b)) and the sum rule is again satisfied but the loss of spectral weight grows 
to ~ 50% at larger couplings (Fig. 2(d)). The breakdown of the one phonon approximation 
is therefore dimension dependent and, by increasing the e-ph coupling, 2D systems seem to 
favour the appearance of multiphononic contributions in the adiabatic regime. 

Figures 3 deal with the intermediate regime u n = J. Electrons are still good excitations 
in the extremely weak coupling and one dimensional case (Fig. 3(a)) with one well defined 
transition at e = —l00.9meV but, at larger couplings (Fig. 3(c)), the sum rule is far from 
being satisfied and a 40% loss of spectral weight is observed together with a strong reduction 
of the main peak height whose position is shifted to e = —119meV. In 2D, the sum rule is 
fulfilled at very weak couplings (Fig. 3(b)) but the appearance of several transitions peaks 
in the range [—110, —90}meV signals the onset of a polaronic state. At larger couplings 
(Fig. 3(d)) the sum rule is strongly violated and the electronic quasiparticle picture totally 
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breaks down. Note that the absorption spectra broaden (both in ID and 2D) by enhancing 
the strength of the e-ph coupling as a consequence of the mixing of electronic states and 
lattice vibrational excitations. Although comparisons with specific data are not possible at 
this stage, broad photoemission spectra are known pj, [T5| , pl}| to be a feature of systems with 
polaronic charge carriers. 

Figures 4 illustrate that electrons are good quasiparticles in the fully antiadiabatic regime 
in one and also in two dimensions. In ID, for both values of the e-ph coupling there is a well 
resolved peak due to the 6— function contribution and located at e = — 98.8meV (Fig. 4(a)) 
and e = —90AmeV (Fig. 4(c)), respectively. The sum rule is well satisfied in Fig. 4(a) 
while a slight loss of spectral weight occurs in Fig. 4(c) being /f^ deA(k = 0, e)/2ir = 0.92. 
In 2D the spectra do not exhibit any relevant change with respect to the corresponding ID 
cases: in Fig. 4(b) a clear transition appears at e = —98.5meV while the peak is located at 
e = —89meV in the moderately weak coupling case of Fig. 4(d). In the antiadiabatic regime 
the main transitions are always due to the first addendum on the r.h.s. of eq.(7). We point 
out that also in the different context of the excitonic spectral function [i5 the disappearance 
of side peaks in the absorption probability had been predicted in the antiadiabatic regime due 
to the fast phonon fluctuations which destroy the high-lying excited states in the potential 
well. 



3. FINAL REMARKS 

The Su-Schrieffer-Heeger tight binding model Hamiltonian has been extended to the 
study of a two dimensional electron-lattice structure. Through a perturbative approach and 
an exact computation of low order diagrams, we study the renormalization of the charge 
carrier effective mass versus the adiabaticity parameter both for a linear chain and for a 
square lattice. In the intermediate regime, where the phonons compete with the electrons 
on the energy scale, we find a sizeable mass enhancement which may be understood as a 
signature of polaron formation. This enhancement is more pronounced in the square lattice 
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mostly at uo w ~ \/2 J. The analysis of the electron spectral function shows that the model 
Hamiltonian hosts quite different behaviors according to the regime set by the adiabatic 
parameter. We have computed the spectral function at the bottom of the band in a number 
of representative cases by varying the strength of the effective coupling and using the spectral 
function sum rule as a testing bench for the reliability of our one phonon approximation. 
While in antiadiabatic conditions the electrons behave as good quasiparticles both in one 
and two dimensions, novel features emerge in the moderately adiabatic and intermediate 
regime where multiphononic terms become appreciable by increasing the strength of the 
e-ph coupling, the spectral weight is progressively spread among several transition peaks 
and the electronic quasiparticle picture is lost. Unlike the Holstein model whose ground 
state polaronic properties are essentially dimension independent, we find that the onset of 
a polaronic state is more likely to occur in 2D than in ID thus confirming the trend of 
the effective mass computation and suggesting that the Su-Schrieffer-Heeger Hamiltonian is 
rather sensitive to dimensionality effects. 
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FIGURES 



FIG. 1. Renormalized masses (in units of bare band electron mass) versus the adiabaticity 
parameter in one and two dimensions, rn^ff 1S ^ ue *° * ne one ph° non self-energy correction. The 
coupling constants are in meV. 

FIG. 2. ID and 2D Electron spectral functions in the adiabatic regime and (a,b) extremely 
weak e-ph coupling; (c,d) moderately weak e-ph coupling. J = O.leV. 

FIG. 3. ID and 2D Electron spectral functions in the intermediate regime and (a,b) extremely 
weak e-ph coupling; (c,d) moderately weak e-ph coupling. J = O.leV. 

FIG. 4. ID and 2D Electron spectral functions in antiadiabatic regime and (a,b) extremely 
weak e-ph coupling; (c,d) moderately weak e-ph coupling. J = O.leV. 
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